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Abstract 

A measurement is presented, using data taken with the HI detector at HERA, of the 
contribution of diffractive interactions to deep-inelastic electron-proton (ep) scatter- 
ing in the kinematic range 8.5 < Q 2 < 50GeV 2 , 2.4 x 10~ 4 < Bjorken-x < 0.0133, 
and 3.7 x 10~ 4 < xjp < 0.043. The diffractive contribution to the proton structure 
function F2(x,Q 2 ) is evaluated as a function of the appropriate deep-inelastic scat- 
tering variables xjp, Q 2 , and j3 (= xjxjp) using a class of deep-inelastic ep scattering 
events with no hadronic energy flow in an interval of pseudo-rapidity adjacent to 
the proton beam direction. The dependence of this contribution on xjp is measured 
to be x^> ra with n = 1.19 ± 0.06(stat.) ± 0.07(syst.) independent of /? and Q 2 , which 
is consistent with both a diffractive interpretation and a factorisable ep diffractive 
cross section. A first measurement of the deep-inelastic structure of the pomeron in 
the form of the Q 2 and (5 dependences of a factorised structure function is presented. 
For all measured /3, this structure function is observed to be consistent with scale 
invar iance. 
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1 Introduction and Formalism 

Recent measurements at HERA of deep-inelastic electron-proton (ep) scattering (DIS) in the 
low Bjorken-x kinematic range 5 < Q 2 < 120 GeV 2 and 10~ 4 < x < 10~ 2 have demonstrated 
the existence of a distinct class of events in which there is no hadronic energy flow in an interval 
of (laboratory frame) pseudo-rapidity rj adjacent to the proton beam direction j|, |j. Our 
present understanding of DIS could be inadequate at low x because additions to the leading 
order QCD-based partonic picture are likely to be substantial. A natural interpretation of these 
so called "rapidity gap" events is based on the hypothesis that the deep-inelastic scattering 
process involves the interaction of the virtual boson probe with a colourless component of the 
proton. Hence there is no chromodynamic radiation in the final state immediately adjacent to 
the direction of the scattered proton or any proton remnant. Observed distributions of such 
events are found to be consistent with simulations based on models in which the virtual boson- 
proton interaction is diffractive ||], that is in which the colourless component of the proton is 
hypothesised to be a pomeron (IP) and the virtual boson-proton interaction may be understood 
as IP exchange. The observation of these rapidity gap events in DIS means that a measurement 
of any short distance sub-structure of this colourless component of the proton is possible, and 
thus, if the process is diffractive, of the IP. An understanding of the sub-structure of this 
colourless component, whether in the form of a partonic interpretation ||, |||, |H| or otherwise, 
is essential for further understanding of the partonic picture of deep-inelastic lepton-nucleon 
scattering. 

In this paper we present a measurement which quantifies the contribution of rapidity gap 
events to the inclusive deep-inelastic structure function F<i of the proton. The results here follow 
our first measurement of this contribution in terms of a "diffractive structure function" F® as 
a function of the two DIS variables x and Q 2 reported in ||. The measurement presented here 
is made as a function of the three kinematic variables (3, Q 2 , and x, or equivalently f3, Q 2 , and 
xjp, which are defined as follows: 



x = 




Xjp = 



q-(P-P' 



) 



Q 2 



= -q 



2 



2q ■ (P - P>) ' 



(1) 



Here q, P and P' are the 4-momenta of the virtual boson, the incident proton, and the final 
state colourless remnant respectively. The latter can be either a nucleon or higher mass baryonic 
excitation, and, if the proton interaction is diffractive, it must be a proton or a proton excitation 
with isospin (1,12,) = (51+3)- Note that 

x = 0x P . (2) 
It is convenient to write xjp and j3 above as 

Q 2 + M 2 x -t Q 2 + M 2 X _ 

X]P ~ Q2 + W 2_ M 2 Q2— -X-XJP/P W 

Q 2 „ Q 2 m 

Q 2 + M\ - t~Q 2 + M\ ~ Xq/IP K ' 

where Mx is the invariant mass of the hadronic system excluding the colourless remnant, t = 
(P — P') 2 is the 4-momentum transfer squared at the incident proton vertex, and W is the 
total hadronic invariant mass. In the kinematic domain of these measurements (M 2 <C Q 2 , 
M 2 <C W 2 ) and if 1 1 | is small (| t |<C Q 2 , \ t |<C M x ) approximating to "the proton's infinite 
momentum frame", xjp may be interpreted as the fraction xjp/ p of the 4-momentum of the 
proton carried by the interacting IP (or meson for non-diffractive contributions), and (3 as the 
fraction x q / jp of the 4-momenta of the IP (or meson) carried by the quark interacting with the 
virtual boson. 

The diffractive structure function F^ 3 \ which is a function of three kinematic variables, is 
derived from a structure function -RP^ which is a function of the four kinematic variables x, Q 2 , 
xp and t. F% is defined by analogy with the decomposition of the unpolarised total ep cross 
section. In the Q 2 range of these measurements, the cross section for the process ep —* eXp, 
where here the final state p specifies both nucleon and higher mass baryon excitation, is assumed 
to be dominated by virtual photon exchange. It is therefore written in terms of two structure 

DU) F u{4) 

functions F 2 and 2 x(i+R D W) m ^ e form (a is here the fine structure constant) 

jW _ 4 ™ 2 L , t. \ F D ^( X O 2 xv t) (5) 

dxdQ 2 dxjpdt ~ xQ* \ V 2[l + R D W(x,Q 2 , XlP ,t)}j 2 [ ^ ' iP ' ' [ 1 

in which y is the usual DIS scaling variable given by y = Q 2 /sx. s is the ep collision centre of 
mass (CM) energy squared. It is convenient to express this cross section in terms of and 
R D ^ because the data available are predominantly at low y so that there is little sensitivity 
to R D ( A \ At fixed xjp and t for the range of y of the results presented here (y < 0.428), 
increases by no more than 17% for < R D ^ < 00. 

The measurement presented here uses data in which the final state colourless remnant is 
not detected. Therefore no accurate determination of t is possible and the measured cross 
section amounts to ^^jq3^j~) from which it is possible to determine only F^ 3 \x, Q 2 , xp) = 

J kP^(x, Q 2 ,xjp, t) dt provided that a particular choice is made for R D ^ and its t dependence. 
The integration is over the range | t m i n \ <\t \<\t \u m where t m i n is a function of Q 2 , W 2 , M x and 
the mass of the colourless remnant, and 1 1 \ u m is specified by the requirement that all particles 
in the colourless remnant remain undetected (see section ^). For this measurement, R D ^ is set 
to for all t and F®^ is evaluated from 



d 3 a e p^ ep x _ 47ra 2 




x,Q 2 ,xip) (6) 



dxdQ 2 dxjp xQ^ \ 

following the original procedure of 0, § . 

The above definition of F^ 3 \x, Q 2 ,xjp) renders quantitative comparison with the structure 
function F%(x, Q 2 ) of the proton straightforward 1|. 



2 HI Apparatus and Kinematic Reconstruction 



The data were taken with the HI detector at the HERA ep collider at DESY in 1993, in which 
26.7 GeV electrons were in head-on collision with 820 GeV protons, yfs = 296 GeV. The HI 



detector is described in more detail in 0, 0, 10, 11]. Here the most important aspects necessary 



to explain the procedures for the analysis are summarised. 

In the following, a coordinate system is used with origin at the interaction point and z axis 
along the proton beam, or forward, direction. The pseudo-rapidity of a final state particle with 
polar angle in the laboratory is then rj = — In tan |. 

Scattered electrons are measured in the backward electromagnetic calorimeter (BEMC, elec- 
tromagnetic energy resolution cr^/E ~ 10%/y/E (GeV) and overall scale known to within 1.7%) 
and in a multi-wire proportional chamber (BPC). 

Charged particles are detected in central and forward tracking detectors (CTD and FTD). 
They consist of drift chambers interspersed with multi-wire proportional chambers for fast 
signals for triggers using charged particles originating from the z range of the event vertex. 

Hadronic energy flow in the event final state is measured in the liquid argon (LAr) calorimeter 
(—1.51 < 7] < 3.65). The hadronic energy resolution is cte/E ~ 50%/y/E (GeV) ©2% (measured 
in a test beam JT2]) and the overall scale is known to within 6%. 



The selection of events in which there is no energy flow adjacent to the direction of any proton 
remnant requires detectors with the best possible coverage in the forward region. The forward 
"plug" calorimeter (PLUG; 3.54 < 77 < 5.08) is used in this analysis to "tag" the production of 
energy above a threshold of 1 GeV. The forward muon detector (FMD) is used in this analysis to 
"tag" the production of hadrons in the forward pseudo-rapidity range 5.0 < r] < 6.6 by detecting 
and reconstructing track segments due to charged particles produced in secondary interactions 
of these hadrons in the collimators, the beam pipe, and adjunct material [13, |ll| ]. 



Two electromagnetic calorimeters (LUMI) situated downstream in the electron beam direc- 
tion measure electrons and photons from the process ep — > epj for the purpose of luminosity 
determination. They are also used in this analysis to help quantify photoproduction background. 

The kinematic quantities x and Q 2 were reconstructed using the "S method" [15, ||, which 
ensures acceptable resolution across the entire kinematic range considered. The procedure, 
which is relatively insensitive to initial state radiative effects, uses information both from the 
scattered electron reconstructed using the BEMC and the BPC and from the hadronic energy 
deposition in the LAr calorimeter and BEMC, together with the event vertex. No attempt was 
made to use energy flow detected in PLUG and FMD for the purposes of reconstructing the 
kinematic variables of events. 

The mass of the hadronic system excluding the colourless remnant, Mx, was determined 
directly from the calorimeter cells in the LAr calorimeter and BEMC and from the event vertex 
determined using the tracks reconstructed in the CTD and in the FTD. (3 was then determined 
using Q 2 from above and Mx substituted in equation @. xjp then follows using equation @. 



3 Data Taking and Event Selection 

The data were obtained with a trigger requiring the presence of a localized energy cluster in the 
BEMC of 4 GeV or more. The selection procedure described in || was used to obtain a sample 
of DIS ep events, which were then constrained to the kinematic ranges 7.5 < Q 2 < 70 GeV 2 and 



Figure 1: Distributions of a) rj ma x measured in the LAr calorimeter for the sample of DIS events; b) the number 
of charged track segments reconstructed in FMD JVfmd per event for DIS events with ri max > 3.2; c) total energy 
deposition per event in the PLUG for DIS events with rj ma x > 3.2; d) hadronic invariant mass of the final state 
system Mx for the selected rapidity gap (RG) events; e) rj ma x for DIS events and for the rapidity gap events 
with the "standard DIS" background subtracted; f) y, g) xjp, and h) (3 for the rapidity gap events with the 
DIS background subtracted; in a), b), c) and d) the expectations of "standard DIS" simulations (CDM, MEPS), 
normalised to the data for rjmax > 3.2, are also shown and the errors are statistical; in e), f), g) and h) the 
expectations of the Monte Carlo simulation (DIFF) used to extract the contribution of rapidity gap events in 
the form of are also shown, and the systematic error due to due to the DIS background subtraction are 

included with the statistical error. 



0.03 < y < 0.7. A total of 16366 events satisfied the selection criteria from an integrated ep 
luminosity of 271 ± 14 nb -1 . 

Figures [j]a), b) and c) show spectra which illustrate the response of HI detectors to these 
DIS events. Superimposed are two simulated expectations, henceforth referred to as "standard 
DIS" , which are based on a partonic picture and which describe the gross features of DIS data 
at HERA g 0]. They are LEPTO |J| (MEPS), which utilises 0(a s ) matrix elements in 



QCD and parton showers, and a combination of LEPTO and ARIADNE |T|] (CDM) in which 
gluons are radiated from a colour dipole. The parton distribution functions for the proton were 
taken from MRS(H) @. The (x,Q 2 ) dependence of F 2 (x,Q 2 ) was adjusted to reproduce the 
most recent measurement of the proton structure function Q. 

For the LAr calorimeter (figure ||a), the distribution of i] max , the pseudo-rapidity of the 
most forward energy cluster |lCj with energy greater than 0.4 GeV, is shown. In all following 
comparisons and analysis, the normalisations of the "standard DIS" simulations are fixed to the 
number of observed events with r] max > 3.2. Comparison with the "standard DIS" expectations 
not only demonstrates the excess of events with an interval of pseudo-rapidity adjacent to the 
proton remnant which is devoid of hadronic energy (rj max < 1.8), as already reported in ||], but 
also illustrates how well the LAr calorimeter response and cluster analysis are understood in 
the forward region (rj max > 1.8). In figure [TJb) the response of the FMD detector, taken to be 
the number of track segments reconstructed in a pair of FMD drift chamber layers (Nfmd), 
for events with rj max > 3.2, is compared with "standard DIS" expectation. For the PLUG 
calorimeter (figure |l|c), the energy spectrum (Eplug) for these same events is shown, together 
with "standard DIS" expectation. It can be seen that the measured responses in HI of FMD 
{Nfmd) and of PLUG (Eplug)-, as well as of the LAr calorimeter (rj max ) in the forward region 
(77 > 3), are all very well reproduced by the "standard DIS" simulations for those DIS events 
with no large rapidity gap. 

A sample of events, henceforth referred to as "rapidity gap" events, in which no final state 
hadronic energy flow is detected adjacent to the proton remnant direction and which is not 
reproduced by the "standard DIS" simulations, was selected by requiring Eplug < lGeV, 
Nfmd < 1> and rj max < 3.2. These threshold values for Eplug and Nfmd were determined us- 
ing randomly triggered events in HI, for which the responses of each are then due predominantly 
to noise coincidences and out of time track segments respectively. In the rapidity gap sample, 
the remaining background of events which can be described by "standard DIS" simulation was 
taken to be that due to the CDM simulation because recent measurements of forward energy 



flow in DIS are known to be well reproduced by it [16]. The uncertainty in this background 
was estimated using the MEPS simulation, which does not describe well the forward energy flow 
in DIS. Crucial to this background subtraction are the efficiencies for rejection of events with 
forward energy flow. They depend on a detailed understanding of the performances of the LAr 
calorimeter, PLUG, and FMD as "taggers" of forward energy flow. A systematic study of the 
response of each detector as a function of the presence or absence of signals in the other two 
and comparison with simulation (CDM and MEPS) yielded a self-consistent set of efficiencies 



with a quantifiable error [13]. 



In figure |l]d) the distributions in M x are shown, for the rapidity gap sample, and for the 
"standard DIS" background (CDM and MEPS). Using Monte Carlo simulation of diffractive DIS 
events (RAPGAP) [§l|, the resolution a{M x )/M x in M x , determined using LAr and BEMC 
calorimeter cells, was found to be < 30% for M x > 6 GeV and to rise to ~ 100% as M x 



decreased from 6GeV to threshold (twice the pion mass). The "standard DIS" background 
remaining after the rapidity gap selection can be seen to vary between 10% and about 50% for 
Mx less than about 25 GeV. In the following quantitative analysis it is subtracted bin by bin, 
and a contribution to the systematic uncertainty in the yield of rapidity gap events is taken 
equal to the difference between the background subtraction calculated using MEPS and that 
calculated using CDM. 

Background, amounting to less than 1% in total, in the rapidity gap sample due to beam 
gas and beam wall interactions, LAr pile up, cosmic rays, two photon processes and Compton 
scattering were considered and removed or corrected for as stated in [|| . The level of background 
due to photoproduction interactions in the rapidity gap sample was determined to be (5 ± 5)% 
in the bin of lowest x and Q 2 , to be less in other bins, and to be in total 1%, both using Monte 



Carlo simulation based on the PYTHIA [22| and POMPYT [23] codes, and using those events 
in the rapidity gap sample in which an electron was detected in the LUMI detectors. 

The rapidity gap sample thus selected and used in the analysis in section ^ below, amounted 
to 1723 events, and 1451 after subtraction of all backgrounds. In figure |l|e) the distribution of 
Vmax is shown both for all DIS events (as in figure |l]a) and for the event sample after subtraction 
of the "standard DIS" background (used to determine - see below). The errors in the 

distribution of rj max in figure |l|e) include the systematic uncertainties for the selection based on 
the forward detectors FMD, PLUG and LAr calorimeter (rj max < 3.2) with the "standard DIS" 
subtraction. 

Below rj max of 1.8 the distribution of the rapidity gap events in r] max is everywhere below 
that for the total DIS sample. This shortfall amounts to 21% in total, of which (7 ± 3)% is due 
to PLUG noise and out of time track segments which fake either or both of energy deposition in 
the PLUG and more than one reconstructed segment in FMD, and (5 ±3)% is due to "standard 
DIS" background. The remaining (9 ± 4)% is attributed to proton dissociation, secondary 
particles from which sometimes generate hits in FMD and PLUG. Monte Carlo simulation of 
DIS diffraction, in which the incident proton dissociates, confirms this interpretation at the 
observed rate if roughly one third of simulated deep-inelastic diffractive events are assigned as 
proton dissociation. 

The use of a forward rapidity gap to select events in which the proton remnant is colourless 
also specifies the kinematic range in t of the measurement presented here. In the selected sample 
of events no particle from the colourless remnant, be it nucleon alone or one of the dissociation 
products of the proton, may be detected in the detector with the most forward sensitivity, 
namely FMD (5.0 < i] < 6.6). It follows that | t | is less than roughly 7 GeV 2 for the range 
x < xjp < 0.05 of this measurement. 



4 Measurement of F^ 3 \/3,Q 2 ,xjp) 



The rapidity gap event sample has been used to measure F®^ 1 (/?, Q 2 , xjp) from the measured 
cross section by means of equation (@). The results were obtained from the data in the form of 
F®^ as a function of /?, Q 2 and x, from which the dependence on j3, Q 2 and xjp follows using 
equation (||). 

Monte Carlo simulations based on deep inelastic electron-pomeron (eiP) scattering (RAP- 
GAP) [21] and on diffractive elastic vector meson electroproduction |24|, |2| (based on |25|, |26|| ) 



were used for the determination of acceptance and experimental bias. Both have been demon- 
strated as necessary to describe well the main features of the data sample Q. Because the 
resulting sample is selected on the basis of laboratory pseudo-rapidity r/ in the forward di- 
rection, it is confined to small values of xjp by virtue of a strong correlation between r\ max 



and xjp |14|| . Therefore is evaluated for the well defined range of the kinematic variable 

xjp < 0.05. 

F®^ was determined using the measured event numbers in three dimensional bins of (3, Q 2 
and x. The bin acceptances were calculated using Monte Carlo simulation to obtain bin averaged 
cross sections fi~3| , 14]. The Monte Carlo simulation modelled deep-inelastic elP scattering, in 
which IP structure was taken to be due to a "hard" quark distribution x q /jp{l — x q /jp) with 
admixtures of a softer gluon distribution (1 — x^/jp) 5 |2l|, [7|, |§. The peripheral t dependence 
for the IP flux "in" the proton followed the parameterisation of 27]. Elastic vector meson 
production was included following [24|. Here Xj/jp, % = q,g are the fractions of the 4-momentum 
of the IP carried by the parton i in the IP's structure. When the extracted bin averaged cross 
sections were included in the Monte Carlo simulation above, the distributions of observables in 
HI were described well. Examples are shown as the curves (DIFF) in figures |]e), f), g) and 
h), where the elP interaction mix is "hard quark" : "soft gluon" : "elastic VMD" 72 : 21 : 7% 
in observed events. Note that the quoted proportion of soft gluon events is for the choice of 
sub-process kinematic cut-off pr > 0.3 GeV. 

The criteria used to choose the kinematic boundaries of each bin included the minimisation 
of effects due to experimental resolution and smearing for (3 and x, and adequate statistics for 
Q 2 . No bin with acceptance less than 30% was included. The application of a selection cut 
Vmax < 3.2 for the sample of rapidity gap events ensured that the invariant mass Mx, and 
therefore (3, was well reconstructed. The diffractive Monte Carlo was used to choose bins in (3 
for which the average error in (3 never exceeded 65% of the width of the bin, and for which the 
systematic shifts from true [3 were always < 5% [13, 14]. None of the bins in (3 included the mass 



range corresponding to the electroproduction of the vector mesons p(770), cj(783) and </>(1020). 
Each bin averaged cross section was then interpolated to a bin centre value in Q 2 ,x) |i~4| ], 
and F®^ was then calculated by straightforward application of equation (^). 

Details concerned with the evaluation of all sources of systematic error in the extraction of 
F®^ are in 13, 14]. They are quoted below as average contributions to F®^ so as to give an 



indication of their relative significance: 



• radiative corrections: 3%; 

• uncertainty in acceptance calculation due to absence of higher order QED terms in Monte 
Carlo simulation: 6%; 

• "standard DIS" background subtraction: 17%; 

• uncertainties in the diffractive simulation for acceptance correction due to a priori igno- 
rance of physics sub-processes: 11% 

• 6% uncertainty in the LAr calorimeter hadronic energy scale: 6%; 

• 1.7% uncertainty in the BEMC electromagnetic energy scale: 4%; 

• 2mrad uncertainty in electron scattering angle 9 e : 4% 

• uncertainty in experimental acceptance due to a priori ignorance of f£* : 8%. 



Having been evaluated for each bin, the above sources of error are combined to form the sys- 
tematic errors quoted in the results in table |l[ To include both proton scattering and proton 
dissociation, an overall correction was made for the shortfall (9 ±4)% attributed to dissociation 
(see section H|). Therefore F®^ is evaluated as the total diffractive contribution to the proton 
structure function F%. 



Figure 2: The diffractive contribution ^ (/3, Q 2 ,xp) to the proton structure function F2 as a function of xp 
for different /3 and Q 2 \ the inner error bar is the statistical error; the full error shows the statistical and systematic 
error added in quadrature; superimposed is the result of the fit establishing a factorisable dependence of the form 
oc x F n (see text). Note that an overall normalisation uncertainty of 8% is not included. 

The results for F®^ 1 are shown as the (/3, Q 2 ,X]p) dependence in figure || and as the (fl, Q 2 ,x) 
dependence in table |l[ Note that the smallest total error on F% is 27%. Therefore there is no 
sensitivity to any choice of R D ^ in equation (|5|) and no significant effect due to the assumption 
R D W = and the use of equation @. An overall uncertainty of 8% is not included in the 
table. It arises from the combination of three uncertainties, in luminosity measurement, in the 
proportion of proton dissociation in rapidity gap events, and in the selection of rapidity gap 
events because of PLUG calorimeter noise and random FMD pairs. The results quoted use the 
"S method" for kinematic reconstruction. The analysis has also been completed using only the 
scattered electron to determine the kinematics and the results agree to well within experimental 
error. Repeating the whole analysis using different cuts, for example no restriction on > 
produced results for F®^ which did not differ to within systematic error from those quoted. 

5 Factorisation of an d Evidence for Diffraction 

Everywhere F^^ is observed to decrease monotonically with increasing xjp in the measured 
range 3 x 10~ 4 < xjp < 0.05. An excellent fit to all data points, irrespective of j3 and Q 2 , is 
obtained assuming a polynomial dependence x~jj? with a single exponent n = 1.19±0.06(sta£.)± 
0.07(syst.), x 2 / d - f - = 32.0/46, 94% C.L. (x 2 /d.f. = 64.5/46, 4% C.L. assuming only statistical 
errors). The observed universal dependence on xp is thus a feature of the measurements at the 
level of statistical accuracy. There is no evidence for any systematic trend in the contributions 
to x 2 as a function of j3 and Q 2 . Using only the reconstructed electron to determine x and Q 2 
an otherwise identical analysis yielded n = 1.30 ± 0.08(stat.) ± 0.16(syst.). 

Such a universal dependence on xjp, independent of f3 and Q 2 , is expected naively if the 
deep-inelastic process responsible for rapidity gap events involves a (colourless) target T in the 
incident proton whose characteristics are not dependent on xjp, and which carries only a small 
fraction of the proton's momentum. The dependence of F^ 3 \p,Q 2 ,xjp) then factorises into 
the product of a universal term (f T / p (xip) oc x F n ), which describes the "flux" of T "in" the 
proton, and a term which describes any structure of T and which is a function of only (3 and 
Q 2 . The electron scattering cross section for the process ep — > eXp, N* etc. can then be written 

where d J^qI x describes the deep-inelastic scattering of the electron off T. 

If the colourless component T of the proton is taken to be well parameterised as a "Reggeised" 
hadronic exchange, then the xp dependence of the flux factor is specified by the leading Regge 
trajectory a(t) in the "asymptotic limit" xjp — > 0, t/s — > 0. The dependence takes the form 
Xjp 2a 1 ' p^, |29|, |3C| ]. In this measurement we perforce integrate over t. To good approximation, 
a{t) may be replaced by at=o because the t distribution of the soft hadronic proton interaction 
is likely to be "peripheral" and the t dependence of a(t) is likely to be slight. 

The observed xjp dependence corresponds to at=o = 1.10 ± Q.03(stat.) ± 0.04(syst.). It is 
therefore consistent with the leading (effective) trajectory which describes phenomenologically 
"soft hadronic" diffractive interactions, namely the IP with a(t) = a*=o + a't and at=o = 1-085, 
a' = 0.25 GeV -2 [31, 32]. If the t distribution of proton diffraction is here like that in soft pp 
collisions, namely peripheral with a dependence e bt with b > 1, then the systematic error, which 
arises by ignoring any t dependence of a(t) of size given by the slope a' = 0.25 GeV -2 of the IP 
trajectory, is < 4%. Note that the leading meson Regge trajectories, the /2(1270)/cj(783) for a 
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Table 1: The diffractive contribution to the proton structure function F2 as a function of j3, Q 2 and x; ST 

is the statistical error, SY is the total systematic error, AC is the (smeared) acceptance, and FB the proportion of 
the signal calculated to be background (see text). Not included in the errors is an overall normalisation uncertainty 
of 8%. 



diffractive-like process with no isospin exchange, and the exchange degenerate p(770)/a2(1320) 
for a process involving isospin exchange, have intercept ctt=o ~ 0.5, giving rise to an xjp de- 
pendence which is much less steep. Parameterisations of pion exchange, either "Reggeised" 
(at=o ~ 0) or otherwise, have only a slight dependence on xjp. Therefore the interpretation of 
the rapidity gap events in DIS at low Bjorken-x as being due in the main either to diffractive 
scattering or to diffractive dissociation of the incident proton can be made without ambiguity. 

The xjp dependence observed here does not specify the M x dependence in the process 
ep — ► epX without a further assumption for the (3 dependence of d/MQ^k-y • However, note 
that X]p = M^ x /s where M£ x is the invariant mass squared recoiling against the proton in the 
inclusive process pe — > p(eX). Therefore the observation of an xp dependence oc x^ 1 in the 
process ep — > epX is equivalent to a dependence on M£ x of (M e 2 x )"" at fixed ep CM energy yfs 
in the inclusive process pe — ► p(eX). It is the distribution in missing invariant mass squared My 
of the form oc (My)"' 2 ^" 1 ', a t =o = 1.085, a' = 0.25 GeV -2 in the soft hadronic interaction 
pp — > pY which is characteristic of high energy inelastic diffraction and which specifies the 
effective IP trajectory a(t) in pp diffractive dissociation []32[| . 



In calculations which use the BFKL QCD formalism of IP dynamics, the IP is found to be 
"hard" with an intercept (at=o = 1 + e, e < 0.5) [33, |34|. Combining in quadrature the statistical 
and systematic errors of the measurement of the xjp dependence above, the intercept cannot 
exceed 1.25 (i.e. e < 0.25) with 99.7% confidence. This does not rule out such a "hard" IP if it 
contributes together with a soft component, such as that described above. A recent calculation, 



Figure 3: Dependence of (f3,Q 2 ) on Q 2 and f3; superimposed on the Q 2 dependence are the results of fits 
at each j3 which assume leading logarithmic scaling violations; the best fit (continuous curve) and the curves 
corresponding to a change of ±1 standard deviation (dashed curves) in the slope are shown; superimposed on 
the (3 dependence are the simplest qq expectation for IP structure - — /?)] (continuous curve), and a constant 
dependence (dashed curve), for which the overall normalisations are determined from fits to the data; also dis- 
played is a dependence of the form [(1 — /3) 5 ] (dotted curve) with arbitrary normalisation. Note that an overall 
normalisation uncertainty of 8% is not included. 

based on a BFKL approach, suggests also that the simplest factorisation with an xjp dependence 
oc ip" may not strictly hold [j35|. The accuracy of this measurement of the xp dependence does 
not exclude such a possibility. 



6 The Deep— Inelastic Structure of the Pomeron 

The observation that rapidity gap events in DIS are dominantly of a diffractive nature, and 
the evidence for a simple factorisation of the structure function F^ 3 \p,Q 2 ,xjp) over values of 
xjp in the range 3 x 1CP 4 < xjp < 0.05 for different values of (3 and Q 2 , lead naturally to the 
interpretation of the (3 and Q 2 dependence of F^ 3 ' as a measure of the deep-inelastic structure 
of diffractive exchange. If such exchange is taken phenomenologically to be described by an 
effective Regge trajectory, this interpretation amounts to the deep-inelastic structure of the IP. 

In figure || are shown the (3 and Q 2 dependences of the integral 



which, assuming factorisation, is proportional to the structure function of the IP. The range 
of integration, xjp L = 3 x 10~ 4 and xjp h = 0.05 is chosen to span the entire xjp range of 
measurements of F^ 3 \ F® is evaluated assuming that the factorisation hypothesis with the 
measured dependence on xp is a good description over this range of integration. The assumption 
is therefore made that, for certain values of (3 and Q 2 , the observed xp dependence is valid 
beyond the measured range in xp, including values which are kinematically inaccessible in this 
measurement. Note therefore that, unlike F®^ or any integration of it over kinematically 
accessible values of (3, Q 2 and xp, F® is not a measure of the diffractive contribution to F2. 
The errors in the results for F% \f3,Q 2 ) include the contributions from the uncertainty above 
(both statistical and systematic) in the "flux" dependence on xp. It is convenient to define 
F-fiP^ 2 ) in this way to avoid the need to specify the normalisation of the unknown diffractive, 
or JP, "flux" in the proton when presenting the data in a form which is suitable for direct 
comparison with theoretical expectation for IP structure. 

There is no evidence for any substantial Q 2 dependence of F% ■ Deep (high Q 2 ) diffractive 
interactions are thus broadly of a scale invariant, and therefore point-like, nature. Furthermore, 
F® is observed to have little dependence on (3 with contributions throughout the measured 
range (0.065 < (3 < 0.65). Therefore the structure which is resolved by the electron in the high 
Q 2 diffractive interaction carries only a fraction of the iP's momentum, implying that the IP 
has its own sub-structure which may presumably be attributed to point-like partons. 

A fit to the hypothesis of a linear log 10 Q 2 dependence demonstrates that the data are 
consistent with zero slope to within ~ 1 standard deviation (<r). The errors on the slopes in 
the fits shown in figure || amount to a = 0.09 (log 10 GeV 2 ) -1 for all (3 except that at f3 = 0.065 
for which a = 0.15(log 10 GeV 2 ) -1 . Therefore the accuracy of this first measurement of the 
Q 2 dependence of F® is such that the possibility of "scaling violations", symptomatic of QCD, 
may be admitted. Furthermore, either they may be (3 dependent and thus characteristic of that 
observed for the structure function of a hadron such as the proton, or they may be insensitive 
to (3 and thus characteristic of the structure function of a "dressed" point-like field quantum 
such as the photon. 




(8) 



Also shown in figure || superimposed on the (5 dependence of F 2 D is the expectation based 
on the simplest qq picture of deep-inelastic diffraction, namely a "hard" quark (3(1 — (5) depen- 
dence [0, [|]. Such a dependence has been justified in an approach in which q, q or both "scatter" 



diffractively off the proton |3^, 35]. In such a picture, diffractive quark scattering at high energy 
amounts to the exchange of two or more gluons and has a leading order energy dependence of 
Regge form [36|. There is a suggestion that the observed (3 dependence may exceed this simple 



qq expectation as — > 0, which is naturally to be expected in any quantum chromodynamic 
interpretation which includes gluons as well as quarks. A fit assuming no dependence of F 2 on 
P is equally acceptable. The ansatz of a soft dependence (1 — /?) 5 is completely ruled out. The 
lack of measurements at large /?, where elastic vector meson production contributes, means that 
no conclusion can be drawn concerning the existence or otherwise of a "super-hard" component 
in JP structure, first suggested by Brandt et al. 



The dependence of F 2 on Q 2 and on P, the appropriate Bjorken-x like variable, resembles 
well established measurements of the structure functions F2 of hadrons in that it is consistent 
with an understanding based on asymptotically free partons. It however contrasts with them 
in that the Bjorken-z dependences of the structure functions F2 of hadrons are all observed 
to decrease substantially with increasing Bjorken-x. This suggests the simplest possible inter- 
pretation for the sub-structure of the IP, namely that due to two "valence-like" partons which 
share the majority of the 4-momentum of the JP, together with modifications at low P due to 
QCD evolution. 



7 Summary and Conclusions 

The contribution to inclusive deep-inelastic electron-proton scattering (DIS) of events, in which 
a region of pseudo-rapidity adjacent to the proton remnant direction is devoid of hadronic energy 
and which is not described in the framework of our present partonic understanding of DIS, has 
been evaluated in the form of a "diffractive structure function" F% (P,Q 2 ,xp). 

The dependence of F 2 on xjp, which may be interpreted as the fraction (xjp/ p ) of the 4— 
momentum carried by the colourless component of the proton with which the electron interacts, 
is measured to be x^T with n = 1.19 ± 0.06(stat.) ± 0. 07 (syst.). This dependence is found to be 
universal, irrespective of the deep-inelastic scattering variables P and Q 2 . This is as expected if 
the electron-proton cross section can be factorised into the product of a term which describes 
the "flux" of colourless component of the proton and a term which corresponds to the cross 
section for the interaction of the latter with the electron. 

In the framework of "Reggeised" hadronic exchange in the t-channel which couples to the 
incident proton, the dependence of F 2 ^ onxp may be interpreted as the intercept at=o = l.lOrb 
0.03(stat.) ±0.04(sysi.) of the leading Regge trajectory. The latter is consistent with that of the 
pomeron trajectory which describes phenomenologically the energy dependence of soft hadronic 
diffraction, and is inconsistent with the intercepts at=o ~ 0.5 of the leading meson Regge 
trajectories. The origin of "rapidity gap events" in deep-inelastic lepton nucleon scattering 
is therefore demonstrated unambiguously to be predominantly diffractive, and the colourless 
component of the proton to be consistent with the pomeron. The precision of the measurement 
of ctt=o does not exclude the possibility of a contribution to deep-inelastic diffractive scattering 
from a harder BFKL motivated pomeron trajectory with a higher intercept. 

The deep-inelastic structure of the pomeron is observed to be consistent with scale invariance 
for all measured values of the appropriate "Bjorken-x" like variable p. F 2 ^ is observed to be 
non-zero for a significant range of values of 0, demonstrating a substantial inelastic contribution. 
The sub-structure which is resolved in diffractive deep-inelastic electron scattering thus carries 



only a fraction of the momentum of the pomeron, meaning, when taken together with the 
observation of broad consistency with scale invariance, that this sub-structure may presumably 
be attributed to point-like partons. 
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